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We have analyzed effects of the hyperfine interaction on electric dipole spin resonance when the 
amplitude of the quantum-dot motion becomes comparable or larger than the quantum dot’s size. 

Away from the well known small-drive regime, the important role played by transverse nuclear fluc¬ 
tuations leads to a gaussian decay with characteristic dependence on drive strength and detuning. A 
characterization of spin-flip gate fidelity, in the presence of such additional drive-dependent dephas¬ 
ing, shows that vanishingly small errors can still be achieved at sufficiently large amplitudes. Based 
on our theory, we analyze recent electric-dipole spin resonance experiments relying on spin-orbit 
interactions or the slanting field of a micromagnet. We find that such experiments are already in a 
regime with significant effects of transverse nuclear fluctuations and the form of decay of the Rabi 
oscillations can be reproduced well by our theory. 

PACS numbers: 85.35.Be,75.75.-c,76.30.-v,03.65.Yz 


Introduction. The interest in coherent manipulation 
of single electron spins has stimulated intense research 
efforts, leading to a great degree of control in a vari¬ 
ety of nanostructures mm- For electrons in quantum 
dots, electron spin resonance (ESR) was first demon¬ 
strated in Ref. 0], However, full electric control of local 
spins might be a better strategy for complex architec¬ 
tures of many quantum dots, envisioned to realize quan¬ 
tum information processing [3j. Thus, electric dipole 
spin resonance (EDSR) was developed relying on either 
spin-orbit couplings urn or the inhomogeneous mag¬ 
netic field induced by a micromagnet [7, [5] . The effec¬ 
tiveness of EDSR is highlighted by recent experiments, 
which could demonstrate Rabi oscillations with frequen¬ 
cies larger than 100 MHz for both approaches u uni- 
To further improve the performance of such spin manip¬ 
ulation schemes, it is important to characterize relevant 
dephasing mechanisms, and especially those which might 
become dominant at strong electric drive. In fact, as it 
will become clear in the following, a sufficiently strong 
drive is able to induce significant and yet unexplored 
modifications on how typical dephasing sources affect 
EDSR. For this reason, while representing the main lim¬ 
itation for accurate spin manipulation, dephasing is still 
poorly understood in large-amplitude regime of EDSR 
(i.e., when the amplitude of motion becomes comparable 
to the quantum dot’s size). 

In this work we will focus on hyperfine interactions, 
which are well known to play an important role in the 
electron spin dynamics of quantum dots. In particular, 
the ESR deplrasing was successfully interpreted in terms 
of a static Overhauser field, with a variance of a few rnT 
in GaAs (3]. The resulting power-law decay and a uni¬ 
versal 7t/4 phase shift of the Rabi oscillations were accu¬ 
rately verified confirming the predominance of nu¬ 
clear spins over other sources of dephasing. While EDSR 


experiments were also generally interpreted assuming a 
power-law decay, the expected t -1 / 2 dependence is vio¬ 
lated at the larger values of the drive mmm • Espe¬ 
cially, Ref. m has demonstrated striking deviations from 
the ESR behavior, including a crossover from power-law 
to gaussian decay. It is also known that the electron mo¬ 
tion, as well as the presence of the drive, can have sub¬ 
stantial effects on spin dynamics and decoherence PJ- 
□21 . These considerations motivate us to provide here a 
detailed characterization of EDSR dephasing induced by 
the hyperfine interaction, paying special attention to the 
large-amplitude regime. As a main objective behind try¬ 
ing to achieve faster Rabi frequencies is to decrease oper¬ 
ation errors, we also establish the limitations on spin-flip 
gate fidelity imposed on EDSR by the hyperfine inter¬ 
action. Finally, we compare our theory with EDSR ex¬ 
periments which, as we will discuss, have very recently 
entered the large-amplitude regime. 

Model. EDSR is induced by a driven periodic dis¬ 
placement of the quantum dot R(t) = e x SR sin cot, which 
we take conventionally along x. For the time-dependent 
wavefunction i/j(r— R(t)), we assume harmonic confine¬ 
ment along the direction of motion (which is applicable 
to both nanowire and lateral quantum dots): 

mr)\ 2 = \ <p(v,z)\ 2 -± r e-^. (1) 
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The spin dynamics can be described with the following 
Hamiltonian: 
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where the first term is the electron Zeeman coupling, 
with e z = gRBB and a the Pauli matrices. The sec¬ 
ond term is the drive, for which we can generally as¬ 
sume b oc SR while other features (e.g., the direction 
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of b) depend on specific details of the spin-orbit cou¬ 
pling or magnetic gradient. The last term in Eq. ([2]) 
is the Fermi contact hyperfine interaction, where n o is 
the nuclear density. R is the spin operator of nucleus 
i, with position r) and coupling A,. The periodic time- 
dependence of Eq. pi) is characterized by Fourier compo¬ 
nents V> m (r) = £ jf /w \ip(f- R(t))\ 2 e~ imut dt, of which 
only the static (to = 0) and resonant (to = ±1) ones 
are of interest here. In fact, in a frame rotating at fre¬ 
quency w ~ e z /h and neglecting fast oscillating terms, 
the transformed Hamiltonian H' reads: 
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where longitudinal/transverse fluctuations are controlled 
by ■i/’o (r) and tpi (t) j respectively. Without loss of gener¬ 
ality, we restrict ourselves to the case b x = b z = 0 HE]: 


H' = 
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where A = e z — Hu is the detuning and h is defined by 
the second line of Eq. ([3]). 

Nuclear fluctuations. On the relatively short time 
scales of the EDSR experiments, it is appropriate to de¬ 
scribe h with a static random classical magnetic field. 
In the lab frame and for infinite-temperature nuclear 
spins, the variance of the Overhauser field is given by 
cr 2 = J2i{^ A i/no) 2 \ip(n)\ A Ii(Ii + l)/3. However, h is for 
a reference frame moving with the dot and rotating at 
frequency w. As a consequence, its statistical proper¬ 
ties differ from the ones in the lab frame. We still have 
(h) = 0, but Eq. |i| implies that (hi), (h1 y ) have an 
interesting dependence on the strength of the drive. For 
finite 5R and |^(r)| 2 as in Eq. |l| we can evaluate (hi), 
(h1 y ) as follows, in terms of Hypergeometric functions: 
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In the above formulas, the only dependence is on 6R/Sx , 
i.e., the amplitude of motion relative to the width of 
the electron wavefunction. Sh z and Sh xy are plotted in 
Fig. 0 a), showing that for SR —> 0 only the longitudi¬ 
nal fluctuations survive. Therefore, in this limit one re¬ 
covers the same behaviour of ESR. At finite 6R/6x, the 
value of Sh z is a decreasing function of amplitude while 
the transverse fluctuations become non-zero and have a 
non-monotonic dependence on 6R/6x. Such transverse 
fluctuations can serve at b = 0 as a driving term PS], 
and in this context were previously discussed through 
an expansion at small SR m or numerical evaluation 




<5R/<5x tir 2 



FIG. 1. (a): Plot of Eqs. (|5| and |6|, which characterize 

the longitudinal/transverse nuclear fluctuations (upper/lower 
curve, respectively), (b): Numerical result for the Rabi oscil¬ 
lations Pflt), averaged over nuclear fluctuations (blue thick 
curves). We used rj = 0.05, A = 0, and SR/5x = 0.1,0.2, ...0.6 
(bottom to top). For clarity, the curves are shifted vertically. 
The red dashed curves are the asymptotic power-law decay 
of Ref. El, valid when SR/Sx < y/2fj = 0.3. The thin grey 
curves are from Eq. 0- (c): Plot of Eq. . Dots are 

decay times from Ref. 1101 , rescaled using TJ — 9 ns and 
rj ~ 0.09. (d): Error in realizing a 7r-rotation, plotted as 

function of 8R/8x at tj = 0.01,0.03,0.05 (bottom to top). 
Numerical results (thick solid) are compared to the asymp¬ 
totic result (thick dashed) and the upper bound (thin red) 
given in Eq. (13 ) [19]. The thin red lines also practically co¬ 
incide with the ESR result, as seen by taking 8h z = 1 and 
8h xy = 0 in Eq. (131. 


m ■ However, both scenarios of EDSR (i.e., based on 
a micromagnet or spin-orbit coupling) are in a physical 
regime distinct from Refs. [T^l - fTK] . because Sh xy is typ¬ 
ically much smaller than the drive b. To see this, we 
notice that Eq. |6j) implies 5h xy < \5R/Sx, thus: 

h x , y <x5h xy ^ SR/Sx __ r] . . 

b~b 2 b/a ~ 2 1 ’ 

In Eq. (0) we defined the useful parameter r\ = aSR/bSx. 
r) is approximately constant (since b oc SR) and is typ¬ 
ically small, according to our later estimates. There¬ 
fore, transverse nuclear fluctuations provide an addi¬ 
tional dephasing mechanism which becomes progressively 
more important, until the maximum effect is reached at 
5R/5x ~ 1.8. We will discuss how the effect of 8h xy 
becomes dominant over Sh z in a regime of sufficiently 
strong EDSR drive, which was already realized by recent 
experiments m- 
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Rabi oscillations. We now use Eqs. (|5j) and © to 
perform a gaussian average Pj,(i) with respect to h of the 
spin-flip probability Pj.(t): 
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Although we cannot provide a general closed-form re¬ 
sult, several relevant features can be explicitly character¬ 
ized. In particular, at sufficiently large drive and detun¬ 
ing we can neglect the components of h perpendicular to 
(b/2)e x + Ae z [see Eq. @1 , to obtain: 
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The Rabi decay time is: 



Tr{ A) 


T* \ ^ I b 2 + 4A 2 

X \)b 2 + ±A 2 {5h z /Sh xy ) 2 ' 
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with TJ = -\/2 h/o the typical inhomogeneous dephasing 
time associated with nuclear spins. Equation ([9]) implies 
a crossover between the ESR power-law decay at weak 
drive to the gaussian decay of the strong-drive regime. 

To exemplify this behavior, we first consider the reso¬ 
nant condition (A = 0) when, besides SR/Sx, the form 
of the decay is determined by T£ and the coefficient g. 
An example of numerical results for P\.{t) is shown in 
Fig- [3b), assuming r) = 0.05. We confirm that the ESR 
power-law decay Gnu is recovered when SR/6x —► 0 
but significant deviations from this known dependence 
appear at larger strength of the drive. In the large-drive 
limit, we find a gaussian decay with no universal 7 t/ 4 
phase shift, as Eq. © becomes an excellent approxima¬ 
tion. This crossover to the strong-drive regime occurs 
when the effect of h x in Eq. Q becomes dominant over 
h z , i.e., bh x h 2 . We estimate the typical values of 
h x , h z using the limit of Sh xy , Sh z at small SR/6x (which 
is justified if g -C 1), i.e., h x ~ \aSR/Sx and h z ~ a. 
This yields the condition: 


SR ^ r— 


( 11 ) 


in good agreement with the numerical results of Fig. 
Of special interest is the decay timescale: 


T R ( a = 0) = 


n 

Sh 


( 12 ) 


xy 


which follows simply from Eq. (10) and is plotted in 
Fig. 0 c). The increasing strength of 5h xy with the drive 
leads to a significantly faster decay, as also seen in the 
time domain results of Fig. ©b). However, it is im¬ 
portant to note that the fidelity of the 7 r rotation grows 
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FIG. 2. Rabi oscillations at finite detuning, (a) and (b) are for 
8R/5x = 0.8 and 0.1, respectively. In both panels, g = 0.05. 
(c) shows the fine cuts of (a) at A/a = 0,2.5,5,... 12.5 
(bottom to top). For clarity, the curves are shifted verti¬ 
cally. The numerical results are virtually indistinguishable 
from Eq. § plotted as grey dashed curves. In panel (d) 
we plot Eq. (10 1 as function of detuning for several values of 
8R/5x = 0.4, 0.8,1.2,1.6 (bottom to top). 


monotonically with SR , as shown in Fig. 0d). For a 
quantitative analysis, we set t = 2nh/b in Eq. Q and 
perform an expansion up to second order in h. After 
statistical averaging, we have: 


1 - max[P;(f)] ~ 


r] 2 {A5h 2 z + TT 2 Sh 2 xl 
{SR/Sx) 2 
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4. 177 2 


{SR/Sx ) 2 ' 


(13) 


Equation (13) includes a contribution proportional to 


Sh xy but the additional dephasing from transverse fluc¬ 
tuations is more than compensated by the decrease of 
Sh 2 and the faster Rabi frequency. Thus, Eq. (13) shows 


that it is always advantageous to apply a stronger drive 
and the effect of the hyperfine interaction on 7r-rotation 
error can be reduced below any desired threshold with a 
sufficiently large SR/Sx [ID] , 

Decay at finite detuning. Considering a finite de¬ 
tuning yields further insight on the role of longitudi¬ 
nal/transverse nuclear fluctuations. The effect of A is 
illustrated in Fig.[2j where panel (c) confirms that Eq. © 
provides an excellent approximation in the strong drive 
regime. A first consequence is that the Rabi oscillations 
approach the “chevron” pattern of Fig. [3a). Further¬ 
more the decay time gets reduced at finite A, which is 
illustrated in Figs. [3c) and (d). 

The dependence of Tr on A has a simple physical ex¬ 
planation, as it can be traced to the difference in strength 
between transverse and longitudinal nuclear fluctuations 
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shown in Fig. [lja). Since Eq. Q implies that a finite 
detuning corresponds to a field along z in the rotating 
frame, the relevant component of the nuclear fluctuations 
(i.e., along the total driving field) becomes a weighted 
average of h x and h z . Since 5h z > Sh xy , the nuclear fluc¬ 
tuations gets enhanced by a finite detuning. As shown 
in Fig. [2])d) , the dependence of Tr on A is particularly 
pronounced at smaller values of SR/Sx. This is natural, 
as the ratio 5h z /5h xy is large in this regime, see Fig.JTJa), 
while the nuclear fluctuations become more isotropic at 
larger SR/5x. Thus, studying the dependence of Tr on A 
allows one to explore how the relative strength of longitu¬ 
dinal/transverse nuclear fluctuations evolves with SR/Sx. 

Stationary limit. We conclude our analysis of the 
Rabi oscillations by commenting briefly on the stationary 
limit P\.{t — > oo), which is a useful quantity to estimate 
a. Common methods rely either on the drive dependence 
at resonance (as done in the ESR experiment of Ref. EH) 
or the linewidth (considering finite detunings). We find 
that, even if the Rabi oscillations are sensitively modified 
by transverse nuclear fluctuations, the effect on P^(oo) is 
negligible in the current experimental regime r) <C 1. A 
significant difference between ESR and EDSR only ap¬ 
pears when 77 ~ 1 (for more detais, see [20| ). Therefore, 
methods to extract a from P^(oo) are still valid for large- 
amplitude EDSR. 

Comparison to experiments. We now discuss the ap¬ 
plication of our theory to available experimental data. In 
Fig. [ 3 ] we show an analysis of the data shown in Fig. 2(c) 
of Ref. [9j, obtained from InSb nanowire dots with strong 
spin-orbit interaction. As seen, our theory is able to re¬ 
produce well the Rabi oscillations and the fit yields val¬ 
ues SR/Sx and b consistent with b oc SR/Sx. For EDSR 
driven by the spin-orbit coupling we have: 


p = 


Iso o' 

2 Sx e z sin 6' 


(14) 


where Iso is the spin-orbit length and 0 the angle be¬ 
tween the spin-orbit field and B. Using Iso = 200 — 300 
nm, Sx ~ 10 — 20 nm [21 j . g = 41, and B = 31.4 rnT 
[D] gives 7) ~ (0.015 — 0.04)/sin 6. Figure [3](b) implies 
77 ~ 0.03, which is consistent with this estimate. 

While Fig. [3jb) has SR/Sx < 0.1, we estimate that 
larger values of SR/Sx were achieved in a recent exper¬ 
iment on GaAs quantum dots EDI- There, the drive is 
based on a micromagnet for which numerical simulations 
give b/SR ~ |.g|/xsx(l mT/nm) [10]. Using the largest 
achieved Rabi frequency / max ~ 120 MHz, |g| ~ 0.4, 
and Sx ~ 35 — 60 nm (corresponding to orbital energies 
~ 0.3— 1 meV), we obtain 5R max /Sx ~ 0.7— 1.2 which is 
relatively close to the condition at which h x fluctuations 
are most effective. On the other hand, the regime of mo¬ 
tional narrowing SR/Sx > 1.8 , 22 l does not appear to be 
within reach of current experiments. For GaAs quantum 
dots a ~ \g\p,B x (1 - 4 mT) [51 [Til [T51 [25] , which allows 
us to estimate a typical range 77 ~ 0.02 — 0 . 1 . 



t (ns) 
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FIG. 3. (a): Fit of our theory to Rabi oscillations of Ref. [9]. 
For each curve, the fit parameters are SR/Sx, a conversion 
factor to current, and b = Snhfu, while <7 = 0.22 ± 0.03 ReV 
is from the experiment [9|. (b): SR/Sx and /r obtained from 
the fits in (a). The error bars are from the uncertainty on o. 
The dashed line is a fit to SR/Sx = C/r. 


Several findings of Ref. m are in good agreement with 
our discussion, including the transition to a “chevron” 
pattern in the strong drive regime | 10 j . As a result, 
our panels (a) and (b) of Fig. [2] are remarkably simi¬ 
lar to the corresponding panels of Ref. ED]. Further¬ 
more, a crossover from power-law decay (for b < \g\figx 
5 mT) to gaussian decay (for b > \g\fiBX 15 mT) was 
observed. The strength of the drive for such crossover 
is compatible with Eq. © , which can be rewritten as 
b > yj2aSx(b/5R) ~ \g\pB x (8—20 mT), using the above 
estimates of a, Sx, and b/SR. We also show in Fig. 0=) 
that Eq. ( 121 is able to reproduce the dependence of Tr 
on the drive strength, with reasonable fit parameters for 
GaAs quantum dots (T| ~ 9 ns, 77 ~ 0.09). 


Conclusion. In conclusion, we have characterized the 
dephasing induced by the hypefine interaction in large- 
amplitude EDSR, and we showed that transverse fluc¬ 
tuations of the Overhauser field are likely to play an 
important role in this regime, recently achieved exper¬ 
imentally. It should be mentioned that also other de¬ 
phasing sources were suggested for EDSR, such as para¬ 
magnetic impurities, charge noise, and photon-assisted 
tunneling [ME]. In the absence of clear evidence (or 
specific predictions) about these alternative mechanisms, 
our theory offers further means to test if nuclear spins are 
indeed the dominant effect, e.g., through a detailed anal¬ 
ysis of Tr as function of both drive and detuning. In fact, 
it is unlikely that other types of dephasing would induce 
the same type of sensitive dependence on A discussed in 
relation to Eq. (10) and Fig. [ 2 ] 


From a more general point of view, our study is help¬ 
ful to assess the limitations to spin manipulation due 
to the quantum-dot motion and the nuclear-spin bath. 
These two aspects are unavoidable for EDSR based on 
III-V semiconductors, in contrast to ESR or spin manip¬ 
ulation based on group-IV materials {24J [55] , Nuclear 
fluctuations, on the other hand, do not represent a fun- 
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damental obstacle to EDSR, since high-fidelity gates can 
be achieved at sufficiently large amplitude. 
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Supplemental material for “Dephasing due to nuclear spins in large-amplitude electric 

dipole spin resonance” 

We discuss here here the behavior of P^(t —> oo), i.e., the stationary value of the spin-flip probability under a 
continuous drive. As a function of detuning, P^(oo) displays a peak around A = 0 and, like in regular ESR (the 77 —0 
limit), a stronger drive leads to a general increase of P^(oo), as well as broadening in A. Representative examples are 
shown if Fig. (Slja). 

We first examine the behavior of the peak value as a function of the drive strength. Taking A = 0 in Eq. (9) of the 
main text gives P^(oo) ~ 1/2 in the large-drive regime, thus the interesting drive dependence of P^(oo) occurs when 
h XiV are smaller than b,h z . By neglecting in Eq. ( 8 ) such transverse nuclear fluctuations, we obtain the following 
expression: 


P;(oo) ~ 


V^b 


A\/2a5h : 


. e f 2 /(8P^) er fc 


2y/2 cr5h z 


(SI) 


which is similar to the ESR result, except that the nuclear fluctuations er are replaced here by a reduced value cr8h z . 
We have checked that Eq. (SI) is in agreement with direct numerical evaluation. As seen in Fig. |Sl[b), deviations 
from the ESR result exist in general but in the current experimental regime can be safely neglected, because they 
only become important when 77 ~ 1 . 

By considering the dependence of P^(oo) on detuning, we can reach a similar conclusions about the EDSR linewidth. 
Figure SI 'c) shows that transverse fluctuations have a significant effect for 77 ~ 1 and b > a, when the linewidth is 
significantly narrower than for ESR. However, very small deviations from the Voigt profile occur for currently more 
typical values 77 <C 1. As stated in the main text, we can thus conclude that P^(oo) is much less sensiteve to the 
physical origin of the drive than the decay of Rabi oscillations. One can see, contrasting Fig. ISlJb) with Fig. 1(b) 
of the main text, that significant deviations from the ESR decay appear at relatively small values of SR/5x 7 when 
P^(oo) has not yet saturated to 1/2. 
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FIG. SI. Characteristics of P\(t —> 00 ). Each family of curves shows results with 77 = 0.05,0.1,0.2,0.5,1. Larger values of 77 
correspond to larger deviations from the 77 = 0.05 results (dashed), which are essentially indistinguishable from 77 = 0. In panel 
(a) we plot Pj_(oo) as function of detuning, at two representative drive strengths b/c r = 0.2 and 2. Panel (b) is the A = 0 value. 
The dots on the dashed curve (77 = 0.05) mark the stationary values for the three lowest curves of Fig. 1(b) of the main text. 
The peak around A = 0 is characterized by the half width at half-maximum plotted in panel (c). 























